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We identify optimal quantum error correction codes for situations that do not admit perfect 
correction. We use entanglement decay as a signature of the correctability of errors, and as specific 
realisations we consider noisy channels with correlated errors on multiple qubits. We present n-qubit 
inequalities that govern the regimes in which different repetition codes are optimal and demonstrate 
significant improvement in performance. 
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I. INTRODUCTION 

The goal of quantum error correction (QEC) is to 
transmit quantum information reliably despite noise 
and decoherence resulting from unavoidable environment 
coupling. On the one hand, QEC is essential for the prac- 
tical implementation of quantum communication proto- 
cols and quantum computational algorithms, but, on the 
other hand, QEC also provides an interesting forum for 
examining the conditions in which generally irreversible 
open system dynamics can be undone. 

Reversibility can only be achieved for specific decoher- 
ence models, and real open system dynamics will typ- 
ically always deviate from such idealized scenarios. In 
practice one may achieve reversibility for the predom- 
inant decoherence effects, but typically not for the re- 
maining subordinate processes. Consequently, if there 
is a clear separation in probability regimes for the var- 
ious detrimental processes, i.e. errors that may occur, 
the best strategy is likely to target perfect correction of 
the more probable errors at the expense of leaving less 
likely errors uncorrected. This is the strategy currently 
employed by most QEC protocols. However, there is a 
plethora of common scenarios in which the simplification 
to separable probability regimes is insufficient, e.g. if, in 
addition to single qubit errors, there are correlated errors 
caused by fluctuations in global control fields, such as the 
laser that creates an optical lattice for neutral atoms [1 
or a field gradient that induces an interaction between N 
ions [2] . Several important steps have already been taken 
in terms of quantifying the impact of correlated noise and 
calculating accuracy thresholds for concatenated codes 
and fault tolerant computation (see for example [3HS])j 
but the simplification remains a persistent problem. In 
the current paper we turn instead to the idea of approx- 
imate QEC: if a clear separation of probability regimes 
is not evident, it might instead be better to abandon the 
goal of perfect correction of the high-probability subset of 
errors [7]. In particular, we use entanglement dynamics 
to determine optimal QEC protocols for the cases where 
a typical separation of probability regimes with dominant 
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single qubit errors is not present. 

Quantum dynamics can be described in terms of quan- 
tum channels f2, i.e. completely positive maps 

N 

g(t) = fi(e(0)) =Y,Ei(t)Q(P)Et(t) , (1) 
f=i 

that connect a system's density matrix g at an initial 
time t = and a later time t. The specificities of the ac- 
tual dynamics are characterized by the N time-dependent 
Kraus operators Ei(t). If the Kraus operators satisfy 

the condition 2i = i^i(0-^iW = 1> then the channel is 
trace preserving and thus deterministic. Conversely, if 
^2i—xEj(t)Ei(t) < 1 then the channel is probabilistic 
and there is a finite probability for a process that is not 
described by the channel to occur. 

In the context of QEC, a deterministic channel is typ- 
ically divided into a correctable and an uncorrectable 
part: 

n t ( g ) = n c (g) + a u ( g ) , (2) 

such that f2 c (f?) has an inverse channel for initial states in 
a well-defined subspace of the full system's Hilbert space. 
This subspace is called the code and f2 c is correctable for 
initial states in the code. If the contributions from Q u 
vanish for t = and result in deviation from trace conser- 
vation of the correctable channel tt c which grows slower 
than linearly with t, then arbitrarily good correction of 
fit can be achieved [8]. 

The question of which states render a given channel 
correctable has a surprisingly simple answer, encapsu- 
lated by the necessary and sufficient conditions for the 
existence of a recovery operation [9j [10]. That is, if and 
only if there is a projector P, i.e. code, such that the 
Kraus operators describing the channel SI satisfy the re- 
lation 

PE\E 3 P = onjP, (3) 

then there is a deterministic quantum channel Q , such 
that — 1 (r2(g») ) = pg with p = TrCl(g) for all states that 
satisfy PgP = g, i.e. all states in the space defined by 
P. Here a is a Hermitian matrix with complex elements 
and the time dependence of the Kraus operators {Ei(t)} 
is left implicit. 
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In the standard implementation of QEC, the condi- 
tions ([3]) are satisfied for the Kraus operators of fl c only, 
while those of fi u reflect the finite (and presumed negligi- 
ble) probability for an uncorrectable error. In the follow- 
ing we consider the scenario that this probability is not 
negligible, and examine the correctability of the full trace 
preserving channel. This more general case occurs when, 
for example, the probability for an error associated with 
the uncorrectable channel grows linearly or faster in t, or 
sufficiently fast repetition of Q,~ l is not possible. 

The task at hand is thus the rigorous identification 
of optimal codes for conditions in which perfect correc- 
tion of the errors is not possible. Such endeavours have 
been attempted previously, and it has been confirmed 
that the performance of quantum codes through the full 
noisy channel can be improved through the relaxation of 
the conditions ^ [7J [TTJ [12] . Approaches to this regime 
of approximate error correction typically require numeri- 
cal optimisation over all possible encodings and recovery 
maps, which is extremely computationally demanding if 
not practically impossible unless one separates the two 
optimisations. The approaches typically use average or 
worst-case fidelity to identify codes with optimal perfor- 
mance. In what follows we shall instead target entangle- 
ment, since its decay is directly connected to reversibil- 
ity, and therefore permits the assessment of code perfor- 
mance without the need to consider recovery operations 



The search for optimal codes thus begins not with the 
projector P, but a set of orthogonal states that span 
the code, i.e. Yli=i \ — P> an d by considering an 

entangled state 



(4) 



between ancillary system with orthonormal states \4n) 
and the system of interest. The system is affected by the 
noisy channel, but there is no dynamics in the ancillary 
system. In other words, Alice keeps the ancilla while 
sending the remaining entangled state through a noisy 
quantum channel to Bob. 

If one has a perfect code, this implies that the state 
of the full composite system can be recovered after the 
effect of the noisy channel through a recovery operation 
that is acting on the system of interest only, but not 
on the ancillary system. In terms of entanglement the- 
ory, this means that the original state can be recovered 
through local operations only. Since the initial state and 
the final state (after the noise and recovery) have the 
same entanglement content (it is the same state), and 
since entanglement cannot increase through local oper- 
ations, this implies that entanglement did not decrease 
through the effect of the noisy channel on the perfect 
code. For imperfect codes therefore, the loss of entan- 
glement due to the noisy channel is a natural tool to 
characterise the correctability of a noisy channel, which 
can now be assessed without the need to consider the re- 



covery operation explicitly. In what follows we present 
a method of targeting entanglement directly to identify 
optimal quantum codes, and demonstrate that such pro- 
cedures provide substantial improvement in performance 
over current QEC implementations. 

The rest of the paper is organised as follows: In section 
[IT]we introduce the regime of approximate quantum error 
correction and additional nomenclature for examples pro- 
vided in the text. We discuss the entanglement dynamics 
in |II A| and demonstrate explicit correlation between the 
rate of entanglement decay and the rate of deviation from 
complete error correction. We describe the code optimi- 



sation process in III and give a simple 3-qubit example 
in |III A| where optimal codes are found to give substan- 
tial improvements. In |IIIB| we present the n-qubit case, 
with an explicit inequality distinguishing the two opti- 
mal regimes. We conclude in IV with some discussion 
and prospects for further work. 



II. APPROXIMATE QUANTUM ERROR 
CORRECTION 



We saw that the possibility of recovering the original 
message by employing a particular quantum code is en- 
capsulated in criterion ([3]) , and we introduced the scenar- 
ios in which no projectors can be found to satisfy those 
conditions. For such cases the deviation from the condi- 
tion of perfect recoverability can be characterized by 



Y,'\r^nXr- with 



2-J 



A, 



PEjEjP - anP, and 
Tr (PEjEjP) 
Tr(P) 



(5) 

(6) 
(7) 



Here the error set {Ei} is now extended to include the 
operators of £l u , thus comprising the full trace-preserving 
channel. If it were possible to correct this extended set 
fully, then Ay — > and ^ would reduce to ^ . Thus as 
an initial hypothesis the condition that S c be minimal for 
a good code seems plausible. This intuitive notion of de- 
viation from (3| was also discussed in [T3] and [HHH], in 
which general and sufficient conditions for approximate 
operator QEC codes were given respectively, along with 
a suggested relationship between the deviation parame- 



ter Aij and the optimal fidelity loss. In section II A we 
demonstrate direct correlation between the deviation and 
the rate of entanglement decay, giving a clear physical in- 
terpretation to the otherwise quite abstract A^. 

For illustrative purposes we shall follow throughout 
the explicit example of a complete dephasing channel, 
through which one typically protects the message - a unit 
of information; a quantum bit - by encoding the message 
into a state with more qubits. The simplest examples of 
such encodings are the archetypal three-qubit repetition 
codes [101 125j which, despite their simplicity, remain im- 
portant test cases for the theoretical development and ex- 
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perimental implementation of quantum computing |18j . 
Error detection with such codes uses the majority-rule 
principle, which fails to produce the correct result if in 
fact two or more errors occurred. Essential to the correc- 
tion procedure for repetition codes is thus the restrictive 
assumption that the channel imparts only single, uncor- 
rected errors on the coded states, an assumption that 
we remove in the following. 

Upon transmission of the message, the noisy channel 
acts on the information according to equation 0. Er- 
rors on qubits are effected by the generators of SU{2) 
(i.e. the Pauli matrices) supplemented with the identity 
operator, and result in bitflip or dephasing errors or a 
combidation of the two. For pure dephasing our building 
blocks consist of the identity / and dephasing error a z : 



I = 



1 
1 



1 

-1 



(8) 



The code space is a tensor product of the individual qubit 
Hilbert spaces, and so the operation elements are con- 
structed accordingly, with the full three-qubit dephasing 
channel given explicitly in Table [I] in the Appendix. We 
assume the probability for single and two-qubit errors are 
independent of the choice of qubits, and since the chan- 
nel is trace preserving p$,pi,p2 and ps respectively give 
the probabilities of errors occurring on zero, one, two or 
three qubits simultaneously. Note that to take account 
of all permutations, the full probability of a single error 
occurring on any one of the three qubits is thus 3pi (or 
3^2 for errors on two of the three qubits) . The procedure 
for constructing the channel for other types and combi- 
nations of errors and codes follows intuitively. 



A. Entanglement dynamics and relationship to 
approximate QEC 

We demonstrate that the rate of entanglement decay 
implies that ^ is a signature of the reversibility of the 
errors effected by the channel by considering the entan- 
gled state Q. The message carried by the system is 
protected by extending the Hilbert space of and for 
three-qubit dephasing the standard code 



p 3 = | + ++)<+ + +| + |-— X- 



with 



|±> = -^(|0>±|1», 



(9) 



(10) 



is one that completely corrects errors on single qubits. 
With this encoding, the entangled state becomes: 



1 



«3 



102,- 



-)), 



(11) 



The corresponding density matrix is p = \ip)(tp\, and the 
entanglement dynamics of p as the coded qubit is trans- 
mitted are calculated using the Lindblad master equation 



(see for example [TU]): 



dp 
dt 



(12) 



Here Lj are the Lindblad operators, which correspond 
to the types of error expected for the channel, and may 
be arbitrary errors in general. For pure dephasing errors 
as in Table |lj the Lindblad operators take the associated 
form L\ = Jyycr z etc, where 71 is the Lindblad 

rate for errors on single qubits, 72 the rate for errors 
on two qubits and so on. Explicit evaluation of (6]) is 
possible after constructing the Kraus operators with ex- 
pressions for the time-dependent probabilities in terms of 
the Lindblad rates (we give an example in the appendix) . 
In turn, the entanglement dynamics can be represented 
using an appropriate measure and in our example we use 
Negativity P) : 



N{p) = \\p T »\\-\ = Y J \W\ 



A; 



(13) 



Here 1 1 . . . 1 1 denotes the trace norm and are the eigenval- 
ues of p Tb which denotes the partial transpose of p over 
subsystem b. 

We confirm that investigating entanglement dynamics 
is a suitable method for assessing the correctability of 
errors via ^ when we examine the relationship between 
([5]) and the decay of negativity of entangled states. At 
t = we start with the maximally entangled state p, for 
which the reduced density matrix p re( i for the code (ob- 
tained by tracing out the ancilla) is a projector. We thus 
explore the relationship between the dynamics that cause 
deterioration of the entanglement of p and the behaviour 
of ([5]) by replacing P in ^ with p red such that 



Pi 



ed,E i Ejp re d — 



^ijPred 



A 



ij ■ 



(14) 



The resulting general expression for Ajj is: 



A; 



ptr Tr(p red ElEjp red ) 

Pred^i^jPred ^7 > Pred- (i-O) 



Tr(p re d) 



By comparing the rate of deviation from ^ (as given 
by (§ using p red in (|l5|)), with the rate of entanglement 
decay for p, we find that the two are strongly correlated, 
as shown in Fig. [T] for the case of a full dephasing channel 
(Table]!]) with three-qubit codes and a particular choice of 
Lindblad rates (71 = 72 = 0.27 c and 73 = j c , where 7 C is 
the arbitrary characteristic rate (s -1 )). The correlation 
is generated by sampling over different codes P' n obtained 
by random unitary transformations of the original code: 



PL = up n u\ 



(16) 



where the U are drawn from a circular unitary ensemble 
[ST] , and also produces the associated new input state 



K) =I®U\n n ). 



(17) 
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The strong correlation that is evident from this sampling 
thus confirms the hypothesis that S c be minimal for good 
codes, since an increase in the rate of change of ^ has 
an associated increase in the rate of entanglement decay. 
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FIG. 1: Scatterplot for 3-qubit dephasing channel shows cor- 
relation between the rate of deviation from |3| as given by Q 
using (15|, versus rate of decay of entanglement of p. Each 



point on the plot represents the value of the initial rate of 
change of the function for a different random unitary trans- 
formation of the code [21]. (1000 points shown). 



III. CODE OPTIMISATION 



The direct correlation demonstrated in Section IH Al 
established that optimal entanglement preservation be- 
tween the ancilla and the system corresponds to the min- 
imum of ([5]), and is the best-case scenario for accurate 
transmission through the quantum channel. We turn 
now to the question of whether we can identify codes 
that produce such optimal behaviour and answer in the 
affirmative. We confirm optimality via numerical opti- 
misation over an initial unitary transformation on the 
code as in (16), i.e. prior to the noisy channel, where 



the parameters are kept entirely general without impos- 
ing the requirement that the code should also completely 
correct single errors, potentially allowing for a tradeoff 
for better performance with the full channel. We begin 
by presenting the details for the explicit 3-qubit dephas- 
ing example, and then present the generalisation to the 
rt-qubit inequalities of the model parameters that deter- 
mine optimal code choice. 



A. Optimal 3-qubit codes for known channels 

Let us examine first the structure and behaviour of 
the standard single-qubit-correcting dephasing code 
The decay behaviour is found by solving the Lindblad 
equation (12) with dephasing rates 71,72 and 73 corre- 



The explicit expressions for the probabilities correspond- 
ing to the three-qubit dephasing channel of Table [I] are 
given in (A-l), and are functions of the Lindblad rates. 



sponding to single, double and triple errors respectively. 



Any particular choice of rates 71, 72, 73 thus governs the 
probabilities and weighting of errors in the channel, and 
also governs the regimes in which different codes are op- 
timal, as will be explained below. 

While the dephasing code being considered is com- 
monly referred to as the code to solve errors on single 
qubits, such a code also satisfies ^ when errors occur 
in other combinations from the full channel (e.g. errors 
occur on zero or two qubits only) since the code cannot 
distinguish between the pairs of Kraus operators that 
contribute to Ay in ( |14[ ). That is, one cannot distin- 
guish between Eq and E7, E% and Eg, E2 and E$, and 
E 3 and £4 of Table [TJ The probabilities corresponding 
to these combinations result in 5 C = 2(^/poP3 + %*JPiPi)- 
This is thus a target that must be beaten for any new 
code to improve on the results of the standard code. 

An exhaustive search of the parameter space of the 
above model shows that choosing a code different from 
the standard ^ can improve performance whenever 
73 > 7i 7 72 1 and that in fact in this regime the decay 
profile always converges to the performance exhibited by 
the code [29j = |0 + +)(0 + +| + |1 - -)(1 - -|. Inter- 
estingly, qubit rotations such as this were also found to 
improve robustness for the noisy evolution of graph states 
in another context [20]. For this new optimal code, the 
Kraus combinations that contribute to violation of ([3]) for 
the full channel are the following elements in Ay: Ej E\, 
e\e^, E\E§ and E^Ej. Since this means there are no 
violating pairs with probabilities P0P2 or pops, such a 
code satisfies Q completely for zero or double, and sep- 
arately zero and triple errors (e.g. if the channel contains 
no errors on single qubits). 

By substituting the corresponding probabilities from 
Table [I] we find S c = 2( v /p Q pi + 2y/pip 2 + ^P2P?) for the 
new code. By comparing expressions for the two codes we 
thus construct an inequality which dictates their regimes 
of optimality: 

2(VpoP3 + 3 VP1P2) > 2( VPoPi + 2 + VP2P3) ■ (18) 

When the inequality is satisfied the code J3 can improve 
upon the original code. From this inequality it is again 
possible to confirm that rates 73 > 71 , 72 yield improve- 
ment of the entanglement decay within the decay lifetime, 
and we provide an example plot of the improvement in 
Figure[2]with the choice of rates 71 = 72 = 0.27 c , 73 = j c 
(for convenience an indication of the corresponding prob- 
ability evolution is given in Table |n] in the appendix) . 
This plot confirms both significant improvement in en- 
tanglement preservation and also demonstrates the op- 
timality of the alternative code which converges on the 
optimised result. The inset is a plot of 5 C , defined in ([5]), 
showing that the steepness of the slope reduces for the 
optimal code, in accordance with the correlation with the 
rate of entanglement decay demonstrated in Figure [l] 
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FIG. 2: Plot shows Negativity ( 13 I versus time for a 3-qubit 
dephasing example with rates 71 = 72 = 0.27 c , 73 = <y c . 
Entanglement decay for the original dephasing code P = \ + 

++} (+ + +| + I } ( 1 shown in purple/below and the 

decay after optimisation over unitary transformation shown in 
orange/above, which is coincident with the decay for the code 

P3 = |0 ++)(0+ + | + |l A significant improvement 

in entanglement preservation is evident. The inset shows the 
corresponding evolution of the violation of the complete error 
correction conditions (J3j) as given by function |5| using (15 1. 
Here the original code is in purple/above and the optimised 
evolution in orange/below. 



Once optimal codes have been identified, the remaining 
step is to implement the appropriate syndrome measure- 
ment and recovery procedure. Consider for example the 
following set of operators: 



A 
A 1 
A 2 

A 3 
with 



I<E) I® I, 

I(g>(T z (g>I, 

{y/2-Zq 2 -q 3 )I®I®a z 

-i(y/3q2+q3-l)(?z ®I®a Zl 



'^,(19) 



< q 2 < §, and < q 3 < 1, 



where the parameters q 2 , q 3 indicate what proportion of 
such errors can be corrected completely. The effect of 
errors from this set can be completely undone with the 
use of the alternative code P3, that is, the conditions ([3| 
are satisfied, which permits the identification of the cor- 
rect recovery procedure in the usual way [10] . Clearly 
the original code cannot satisfy ^ for this choice since 
it includes both two- and three-qubit errors along with 
the identity. We can see that if q 3 = 1, we can correct 
simultaneous three-qubit errors completely, while sacri- 
ficing correction on simultaneous errors on qubits two 
and three, with a similar relationship governing A 2 . For 
the choice of operators ( 19 ) it is possible to find the set 



ensure maximal recoverability of errors in the channel: 
P2 > P3 => <72 = § , 93 = 
P2 < P3 <72 = l,q3 = 1- 



The choice ( 19 1 is thus an example where errors on multi- 



ple qubits can be corrected completely by the alternative 
code and not the original. 

The above method of identifying optimal codes also 
holds for the bitflip channel, which is unitarily equivalent 
to the dephasing channel. For bitflips the Lindblad equa- 
tion (12 1 retains terms with the identity, which implies 



that the evolution of the Kraus operators will depend on 
all four Lindblad rates 70, 71, 72 and 73. Similar improve- 
ment to the dephasing case can thus be found for bitflips, 
with the same relationship between the rates separating 
the two regimes of optimal performance: 73 > 71, 72, i.e. 
the division is independent of 70, with the original code 
being (|000)(000| + | 111) (111 | ) , and alternative optimal 
code (| +00)(+00| + I - 11><— 11|). 



B. Optimal n-qubit codes for known channels 

The generalisation to n-qubit repetition codes proceeds 
as above with appropriate extensions of the error chan- 
nel and associated Lindblad rates. While in general it 
is desirable to keep the number of qubits to a minimum, 
analysis of n-qubit coded states gives further insight into 
the structures leading to optimality. Such higher-qubit 
systems are frequently needed in practice, and the analy- 
sis below shows that substantial improvement of the en- 
tanglement decay profile is also possible in these cases. 

As in the 3-qubit case above, optimisation over the pa- 
rameters of an initial unitary transformation of extended 
qubit systems reveals convergence to the same decay pro- 
file as given by two classes of optimal codes: the original 
repetition code 



P n = I +® n ){+® n I 
and the new code 



K = |o, 



<gm— 1 1 



0n- 



(20) 



- 1 )(l,- 55n - 1 |.(21) 



This latter notation indicates the coded state contains a 
single |0) in the logical code for |+), supplemented with 
(71 — 1) qubits in |+), and the corresponding arrangement 
of a single |1) and (n— 1) qubits in |— ) in the logical code 
for |— ). When the probabilities of the different classes of 
errors (single, double etc.) are independent of the choice 
of qubits as is the case in our channels (e.g. Table |l| , 
it does not matter which qubit of the code contains the 
rotation. 



Similarly to equation (18), we can derive an n-qubit 



inequality which separates the regimes in which the two 
n-qubit codes (20) and (21) produce optimal (slowest) 
rate of entanglement decay: 



of conditions on the choice of q 2 and 93 that will always 



i=0 



71-1 



VPiP 



(22) 



() 



where (™) are binomial coefficients. When the inequality 
is satisfied, the new code (21) improves upon the perfor- 



mance of ( 20 ) . 

The utility of the method is underlined by consider- 
ing an ?i = 4 example in Fig. |3j in which we compare 
the 4-qubit code (20) correcting single qubit errors, the 

-®3 1 _,_ I ■ C 



_®3 



>( + ,- 



-)(- 



which 



code P 4 = | 

corrects errors on four qubits, and the code (21) with 
4 qubits, which also corrects correlated errors on all 
four, but was found in the text to be optimal also for 
the intermediate range. We have chosen Lindblad rates 
7i = 0.27c, 72 = 0.37c, 73 = 0.±7 G , 74 = 27 c (associated 
indicative probability evolution is given in Table III in the 



appendix), and it is evident that the new code (21) sig- 



nificantly outperforms both of the previous approaches. 
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FIG. 3: Plot shows Negativity (131 versus time for a 
4-qubit example with rates 71 = 0.27 c , 72 = 0.37 c , 
73 = 0.17c 74 = 27 c . An indication of the evolu- 
tion of the probabilities of the associated errors is given 
in Table |III| Entanglement decay for initial codes Pa 
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shown in purple/below, code 



®3 



>(+, 



+ + 



g>3 



->(+ 



®3 



— j in blue/middle, 

®3\/i _®3| 



P* = l+, 

and code P' A = |0, +® 3 )<0, +^\ + |1, -® 3 )<1, -~~| m or- 
ange/above. A significant improvement in entanglement 
preservation is evident. The inset shows the corresponding 
evolution of the violation of the complete error correction con- 
ditions (|3| as given by function |5| using ( 15 \ with purple as 
the top line, blue in the middle and orange below. 



profiles of different codes thus enables performance evalu- 
ation of different codes and gives a clear indication of how 
correctable errors of a particular channel are, without the 
need to consider recovery procedures directly. Such tools 
are invaluable in determining the choice of code in prac- 
tice, where every additional qubit required for quantum 
coding is a significant barrier to implementation. 

We have presented the explicit details for optimal per- 
formance in cases where there is a single type of error, 
such as dephasing. The same approach, however, also 
applies directly to situations in which different types of 
errors affect the qubits. The presence of multiple types 
of errors (i.e. bitflip in addition to dephasing) typically 
requires additional overhead in the number of qubits that 
realise a code. Whereas the minimal number of qubits 
to define a code to protect against arbitrary single-qubit 
errors is 5 |22 , 23j , practical realizations often use even 
more qubits such as the 7-qubit Steane code [21] or the 
9-qubit Shor code [25]. Since rates for bit-flips are typ- 
ically substantially lower than those for dephasing, ap- 
proximate error correction certainly defines pathways to 
work with fewer qubits than necessary for perfect correc- 
tion while improving overall performance. 

For general quantum information processing, it is es- 
sential to define tools that work for arbitrary states. In 
more specific applications such as quantum simulations 
[2"B"H2"5] . it might be desirable to stress the correctability 
of certain states more than others. Despite the tremen- 



dous usefulness of inequality (22), it has no generaliza- 
tion to cases where the projector P in ^ is replaced by a 
density matrix g with different weights for different com- 
ponents, which would permit preference to be given to 
certain states. In that case, optimization of reversibility 
via the violation of (T31) is no longer possible, but defining 
an entangled state \k) between code space and ancilla still 
enables the assessment of reversibility via entanglement 
decay. 
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IV. OUTLOOK 
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Appendix A: Supplementary information for 
dephasing examples 



The operation elements of the full dephasing channel 
are given in Table [Tj 
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TABLE I: Table of error set for full dephasing channel for 
three-qubit codes. 



Explicit expressions for the time-dependent probabili- 
ties of errors occurring in noisy channels can be found by 
equating the Lindblad form with the Bloch representa- 
tion of p in the standard way (see e.g. [TO])- For the full 
dephasing channel given in Table [TJ the explicit expres- 
sions for the probabilities in terms of the Lindblad rates 
are: 
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Indicative values of the probability evolution of (A-l I for 
a particular choice of Lindblad rates where improvements 
can be made to the standard dephasing code are given in 
Table HI 
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TABLE II: Indicative approximate values of the probability 



evolution of (A-l I for a 3-qubit repetition code sent through 
a full dephasing channel (Tab le [I] ) with rates 71 = 72 = 0.27 c , 
73 = 7c as represented in Fig. |2( Here the probabilities of sin- 
gle, double and triple errors are approximately equal for short 
times, and single and double errors begin to dominate at long 
times although the contribution from 713 remains significant. 



r 



Indicative values of the probability evolution for n = 4 
with Lindblad rates 71 = 0.2j Cl 72 = 0.37 c , 73 = 0.l7 c , 



74 = 27 c are given in Table III 
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TABLE III: Indicative approximate values of the probability 
evolution according to a 4-qubit example with rates 71 = 
0.27c, 72 = 0.37c, 73 = 0.l7c 74 = 27c corresponding to 
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